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ABSTRACT 

Context. We present the discovery of two transiting extrasolar planets by the satellite CoRoT. 

Aims. We aim at a characterization of the planetary bulk parameters, which allow us to further investigate the formation and evolution 
of the planetary systems and the main properties of the host stars. 

Methods. We used the transit light curve to characterize the planetary parameters relative to the stellar parameters. The analysis of 
HARPS spectra established the planetary nature of the detections, providing their masses. Further photometric and spectroscopic 
ground-based observations provided stellar parameters (log g,T^ff,v sin i) to characterize the host stars. Our model takes the geometry 
of the transit to constrain the stellar density into account, which when linked to stellar evolutionary models, determines the bulk 
parameters of the star. Because of the asymmetric shape of the light curve of one of the planets, we had to include the possibility in 
our model that the stellar surface was not strictly spherical. 

Results. We present the planetary parameters of CoRoT-28b, a Jupiter-sized planet (mass 0.484±0.087Mjup; radius 0.955 ±0.066Rjup) 
orbiting an evolved star with an orbital period of 5.208 66 ± 0.000 34 days, and CoRoT-29b, another Jupiter-sized planet (mass 
0.85 ± 0.20Mjup; radius 0.90 ± 0.16Rjup) orbiting an oblate star with an orbital period of 2.850 548 7 ± 0.000004 1 days. The reason 
behind the asymmetry of the transit shape is not understood at this point. 

Conclusions. These two new planetary systems have very interesting properties and deserve further study, particularly in the case of 
the star CoRoT-29. 

Key words, stars: planetary systems - techniques: photometry - techniques: radial velocities - techniques: spectroscopic 


1. Introduction 


Spacebo rne surveys of transiting extrasolar planets like 
CoRoT dBaglin et al.ll2006l) or Kepler dBorucki et alJl2010h have 
provided crucial evidence of the interactions between stars and 
planets and their common evolution. In this paper we report the 
discovery and characterization of two systems with interesting 
properties. 


* The CoRoT space mission, launched on December 27th 2006, was 
developed and is operated by CNES, with the contribution of Austria, 
Belgium, Brazil, ESA (RSSD and Science Programme), Germany, and 
Spain. Based on observations obtained with the Nordic Optical Tele¬ 
scope, operated on the island of La Palma jointly by Denmark, Pin- 
land, Iceland, Norway, and Sweden, in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astrofisica de Canarias, in 
time allocated by OPTICON and the Spanish Time Allocation Com¬ 
mittee (CAT). The research leading to these results has received fund¬ 
ing from the European Community’s Seventh Framework Programme 
(FP7/2007-2013) under grant agreement number RG226604 (OPTI¬ 
CON). This work makes use of observations from the LCOGT network. 


CoRoT-28b is a hot Jupiter orbiting an evolved star. We in¬ 
vestigated the stellar and planetary properties, in particular the 
rotation and metallicity of the star, and found a lithium abun¬ 
dance higher than expected for its evolutionary state. 


CoRoT-29b is a bit of a riddle. The transit discovered by 
CoRoT is significantly asymmetric, a signature that the host star 
could have a non-spherical shape. The asymmetry measured 
by CoRoT has been confirmed by ground-based observations 
at several observatories. We developed a code that successfully 
models the transit light curve assuming that the star is oblated. 
The deformation of the star implied by our model, however, is 
far too large compared to our current understanding of stellar 
interiors. Other possible scenarios are bands of stellar spots on 
the surface of the star, which if properly placed, could also repro¬ 
duce the transit light curve measured. We show the observational 
evidence and our modelling efforts, but we leave open the ques¬ 
tion of the origin of the asymmetry. 
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Table 1. IDs, coordinates, and magnitudes of CoRoT-28 and CoRoT-29. 


CoRoT name 
window ID 
CoRoT ID 
USNO-A2 ID 
USNO-Bl ID 
2MASS ID 

Coordinates 

CoRoT-28b 
LRc08 L2 0275 
652180991 
0900-13187173 
0955-0374996 
18344520-H0534254 

CoRoT-29b 
LRc08 L2 3502 
630831435 
0900-13220969 
0964-0373597 
18353650+0628467 

RA (J2000) 

18h 34m 45.21s 

18h 35m 36.50s 

Dec (J2000) 

5°34'25.72" 

6°28'46.99" 

Magnitudes 

Lilter 

values 

B“ 

14.994 + 0.033 

16.704 + 0.105 

Y« 

13.775 + 0.022 

15.560 + 0.069 


13.246 + 0.012 

15.090 + 0.060 

pa 

12.627 + 0.036 

14.515 + 0.025 


11.601 + 0.023 

13.565 + 0.021 


11.143 + 0.025 

13.120 + 0.021 


11.032 + 0.025 

13.048 + 0.031 

WR 

10.912 + 0.022 

12.974 + 0.030 

W2'^ 

10.987 + 0.022 

13.064 + 0.036 

W3'^ 

11.866 + 0.314 

12.088 + 0.295 

W4^ 

9.031 (upper limit) 

9.031 (upper limit) 


Notes. Provided by ExoDat <Deleuil et alj|20^: from 2MASS 
catalogue; WISE photometry from !Wright et alJ OOlOli . 


CoRoT-28b - CoRoT raw light curve 
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Fig. 2. Raw light curve of CoRoT-28. The original data points, sampled 
at 32s, are shown in blue. We have binned the data to 512s bins to guide 
the eye. 

ing the planetary parameters from the light-curve fit. Standard 
methods consist in the subtraction of a low-level polynomial (in 
our case, a parabola) to an interval around each transit with a 
length of a few transit durations (in our case, three). This method 
was applied successfully with the light curve of CoRoT-28, and 
we then proceeded to the modelling of the planetary parameters 
as described in subsection l4.ll 


2. CoRoT observations and data reduction 


2.1. General description 

The satellite CoRoT observed the field LRc08 (pointing coordi¬ 
nates 18h 28m 34.7s, 5°36'0.0") in 2011 between 6 July and 30 
September. The CoRoT data for the run LRc08 has been public 
since 7 June 2013 and can be obtained through the IAS CoRoT 
Archiv43- 

The coordinates, identification labels, and magnitudes of the 
stars CoRoT-28 and CoRoT-29 are given in Table [1] CoRoT- 
28 is a relatively bright target and it was assigned an observing 
cadence rate of 32s from the start of the observations, collect¬ 
ing 225 087 measurements with a chromatic mask. CoRoT-29 
is a fainter target and was assigned the cadence rate of 512s. 
When the first transits were discovered by the Alarm Mode, 
the target was assigned a sampling rate of 32s, collecting in 
total 121 288 measurements with a monochromatic mask. The 
region around both targets, indicating the position of the near¬ 
est contaminants and the orientation of the masks, is shown 
in Fig. [T] F or an overview of the CoRoT observing modes , 
please refer to ^isnard & Auvergn3(l2006h ; |Barge et al.l(l2008l ): 
I Auvergne et akri 20091). 


2.2. CoRoT-28b 

The raw light curve of CoRoT-28 is shown in Fig. |2] We treated 
the EN2 leveQ light curve of CoRoT-28 observed by CoRoT in 
its version 3.0. The removal of the influence of instrumental sys- 
tematics and stellar variability is a necessary step before obtain- 


* http://idoc-corot.ias.u-psud.fr/ 

^ For more inform ation about the CoRoT data products, please see 
iBaudin et alJ l2006l) . 


2.3. CoRoT-29b 

In the case of CoRoT-29, the standard reduction procedures de¬ 
scribed in Section l2~2l revealed a feature with the shape of non¬ 
zero slope in the bottom of the transit shape, which is otherwise 
expected to be reasonably flat. The slope of the bottom part of 
the transit is non-zero with 95% confidence level in the 32s sam¬ 
pled data by CoRoT, and only non-zero with Icr significance in 
the 512s sampled data set. But at the same time, CoRoT-29 is a 
fainter star than CoRoT-28 and its light curve has more noise at 
the timescales comparable with the transit length. We therefore 
decided to use a more robust filtering for this light curve. _ 

We followed the procedure described in lAlanini & Aigraini 
(l2009ll . which iteratively fits the planetary parameters, subtracts 
the planetary model from the data, creates a model of the stel¬ 
lar variability, then subtracts the stellar variability model from 
the original data, and then fits the planetary parameters again 
until the method converges to a final set of planetary parame¬ 
ters and stellar activity pattern. In the first iteration, we used a 
sliding time window of 0.5 days in length to deal with the varia¬ 
tions due to activity. For the interpolation, we removed the points 
measured during the transit. But for the next iteration steps we 
divided the light curve in intervals of 0.5 days and fitted each in¬ 
terval with a Legendre polynomials of order 5. This was chosen 
because a sliding window with short timescales might have an 
influence on sharp features in the light curve, such as the transit 
ingress or egress, and the method migh t fail to converge (see also 
the discussion in iBonomo et al.ll2012h . 

We tried with intervals from 0.5 to 1.5 days in steps of 0.1 
days and though the removal of the systematics is not equally 
satisfactory for each interval, the non-zero slope in the bottom 
transit of CoRoT-29b remained. A Legendre polynomial of or¬ 
der 5 within a 0.5 day window provides a smooth fit to the data 
and it is able to remove the contributions from the out-of-transit, 
spot-induced stellar activity and low-frequency instrumental sys- 
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2'x2' POSS image centered on Corot-28 
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2'x2' POSS image centered on Corot-29 



#: Rmag 
0: 15.2 
1: 14.6 
2: 16.7 
3: 17.3 
4: 17.1 
5: 16.4 
6: 15.3 
7: 15.9 
8: 17.8 
9: 17.9 
10: 17.3 
11: 18.3 
12: 17.9 
13: 15.8 
14: 16.0 
15: 16.8 
16: 14.7 
17: 16.8 
18: 15.9 
19: 17.9 
20: 14.6 


Fig. 1. Second Palomar Observatory Sky Survey (POSS II. lReid et alJll99lh image around the CoRoT-28 (left) and CoRoT-29 (right) targets. We 
have over-plotted the CoRoT masks used in the observations and we give identification labels for the brightest, closest background stars. The mask 
of CoRoT-28 shows the subdivision in the three submasks, named from north to south; red, green, and blue. 


tematics present in the data. The regions of the light curve sam¬ 
pled at 512s and at 32s were treated as independent data sets. 
Our method converged after three iterations (one with the slid¬ 
ing window and two with the Legendre polynomials, see Fig.[3]l. 


2.4. Background correction of CoRoT-29b 

Comparing the transits of the section of the CoRoT light curve 
sampled at 512s and at 32s, we were able to identify a change 
in the o bserved transit depth, which required further study (see 
Fig. |4]i. lAuvergne et al.l (l2009l) described the photometric capa¬ 
bilities of the CoRoT satellite in detail. However some of those 
capabilities have slightly degraded with the ageing of the instru¬ 
ment. One effect is the increase of the dark current of the CCD, 
which has an important impact on the on-board photometry. The 
dark current of the CCD was modelled as a uniform value cal¬ 
culated in 196 (14x14) windows distributed along the CCD; 3/4 
of them sampled at 512s and 1/4 of them sampled at 32s. The 
correction applied to the 32s and 512s light curves was the me¬ 
dian value of the 32s and 512s background windows, respec¬ 
tively, to mitigate the impact of hot pixels and cosmic ray im¬ 
pacts. However, with the ageing of the CCD, the uniform model 
for the dark current is no longer a faithful description of the phe¬ 
nomenon. We observed a gradient along the Y axis of the CCD 
whose slope has increased with time. Moreover, the distribution 
of background windows sampled at 512s and 32s is not uniform, 
which introduces an additional bias in the background correc¬ 
tion. The median background level at 32s sampling is calculated 
in a region of the CCD that has in average a lower Y value than 
the average of the windows sampled at 512s, which introduces a 
visible difference in the mean level of the light curve at the two 
different sampling rates. This effect is common to all light curves 
on the CCD during this particular run, but it cannot be seen in 
the light curve of CoRoT-28, as this target was only sampled at 
the 32s cadence. This effect can be mitigated by modifying the 
corresponding background correction for each star as a function 
of its position on the CCD, a feature that will be introduced in 
future releases of CoRoT data. 

The CoRoT pipeline is currently implementing an improved 
background correction to mitigate this effect. We used the data 


of this improved pipeline (version 3.^. However, the final plan¬ 
etary parameters determined when using the previous (version 
3.0) or the new data (version 3.5) are not significantly different. 
Our modelling fits the background contamination as a free pa¬ 
rameter, and can deal with the residuals left by the background 
correction from the pipeline. The incomplete background correc¬ 
tion pipeline dilutes the transit signal as if there was an additional 
contribution of light from a background star, a well-known effect 
that can be corrected with standard procedures. 

3. Ground-based observations 

The main objective of ground-based photometric follow-up is 
to check whether the observed transit features occur on the 
targe t star or might a rise on a nearby eclipsing-binary sys¬ 
tem dPeeg et al.ll2009l) . Then, ideally, low-resolution spectro¬ 
scopic measurements characterize the nature of the stellar tar¬ 
get before high-resolution radial velocity measurements inde¬ 
pendently confirm the planetary nature of the target providing 
a value of its mass (iMoutou et al.ll200^ . This section describes 
the ground-based observations carried out to characterize these 
CoRoT targets. 

3.1. Photometric measurements 
CoRoT-28b 

Photometric follow-up of CoRoT-28b was performed with 
the IAC80 (Tenerife) and the Euler (La Silla) telescopes. Pho¬ 
tometric data were acquired at the 80cm telescope of the lAC 
on Tenerife (IAC80). Measurements during the transit (on) ob¬ 
tained on 9 October 2011 and out-of-transit (off) data taken on 
25 October 2011 showed an on-off brightness variation of 0.3%, 
but at a low level of confidence. Each observing run lasted 33 
minutes. Conversely, relevant brightness variations could be ex¬ 
cluded with high confidence on the neighboring stars (stars 2, 
3, 6, 13, and 16 in Fig. [TJ. A star some 17"SSW of the target, 
on which photometry was performed, was found to be 5.8 mag 

^ The version 3.0 of the data is already publicly accessible through the 
IAS CoRoT Archive (http://idoc-corot.ias.u-psud.fr/). The version 3.5 
of the data will be publicly available in coming months. 
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Fig. 4. Change in the transit depth of CoRoT-29b observed in the light 
curves sampled at 512s and 32s (version 3.0 of the data). 
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on 29 August (on) and on 2 September (off), no eclipses were 
found, neither on-target nor on any nearby stars. Given that the 
uncertainty in the ephemeris at the time of both observations was 
only 4 minutes, we can therefore conclude that the transits must 
occur on the target. 

CoRoT-29b 
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Fig. 3. Iterative filtering applied to the light curve of CoRoT-29b, con¬ 
verging in three steps (labelled A, B, and C). Top of each panel shows 
the light curve (normalized flux vs CoRoT Julian Date, which is HJD 
- 2451545.0), the bottom figures on each panel show the folded light 
curves (normalized flux vs orbital phase). 


fainter than the target and hence too faint to cause a false alarm; 
the recognizable stars that were closest to the target were even 
fainter (about 4"N and 5"S; no photometry could be performed 
on them). In the data from the Euler telescope, taken in 2011 


Photometric follow-up observations of CoRoT-29b were per¬ 
formed on several occasions. A full transit was obtained on 28 
May 2012 with the 2m Faulkes Telescope North (FTN). The 
analysis of this data showed a clear transit on CoRoT-29 itself. 
There are several nearby stars which fall within the same PSF 
and aperture in the CoRoT data, contaminating flux of the mask 
of CoRoT-29. Most of this contamination is due to a similarly 
bright star at a distance of 10"NW which can be well separated 
in the ground based images. Therefore the real depth of the tran¬ 
sit, as it is observed from ground is 1.5% (measured from FTN 
observations) instead of the 0.6% in the CoRoT light curve, be¬ 
cause the contaminating star accounts for about 50% of the flux 
in the CoRoT mask (see Table 0. The FTN observations show 
also the non-flat feature at the bottom of the transit observed by 
CoRoT (see Fig.|5]l. Given the on-target detection of the transit, 
contamination from eclipsing binaries at distances larger than 
2" could be excluded as a source of false alarm. 

Further ON-OFF photometry undertaken by the 3.5m CFHT 
on 10 May 2013 also showed an on-target transit with a depth of 
1.4+ 0.3%. 

A full transit of CoRoT-29b was observed with IAC80 in 
the night of 14 July 2014. The transit centre time, at BID 
2456,853.435 ± 0.005, was used to refine the ephemeris in Ta¬ 
ble |6] It also shows some degree of asymmetry, though it is dif¬ 
ficult to assess it with the same confidence level as for FTN data 
(see Fig. 01. The ingress is partially missing, complicating the 
assessment of any possible out-of-transit slope, either caused by 
stellar activity or by residuals of the data reduction (airmass cor¬ 
rection, instrumentals, etc.). 

3.2. Radial velocity measurements 
CoRoT-28b 

Radial velocity measuremen t s were first obtained with the 
SOPHIE (iPerruchot et al.l 120081 : iBouchv et al.l l2009l) spectro¬ 
graph on the 193 cm telescope at the Observatoire de Haute 
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Fig. 5. Superposition of the transits of CoRoT-29b observed by 
CoRoT (small gray points) in 2011 and the transit observed by FTN 
in 2012 (large blue points). This is not a fit. In this Figure we just su¬ 
perpose the FTN and the CoRoT data corrected from the contamination 
measured in the light curve. The contamination value was obtained ex¬ 
clusively from the fit to the CoRoT data. See text for discussion. 
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Fig. 6. Superposition of the transits of CoRoT-29b observed by 
CoRoT (small gray points) in 2011 and the transit observed by IAC80 
in 2014 (large blue points). Same conditions as in Fig.[5] 


Provence (France) and showed variations in phase with the 
CoRoT ephemeris, establishing the planetary nature of the can¬ 
didate. A total of 25 radial velocity measurements were acquired 
in High Efficiency mode (R - 40 000) with SOPHIE between 
27 August 2011 and 17 August 2012. An additional seven mea¬ 
surements between 7 July 2012 and 10 August 2013 were per¬ 
formed using t he HARPS spectrop aph in standard HAM mode 
{R = 115000) dMavor et al.ll2003lb mounted on the 3.6 m ESO 
telescope at La Silla Observatory (Chile) as part of the ESO 
large program 188.C-0779. Einally, ten spectra of CoRoT-28 
were taken at different epochs between June and July 2012 using 
the Elbre-fed fchelle Spe ctrograph (EIESt l^andsen & Lindber^ 
ll999l:Pfilting et alJl2014l) mounted at the 2.56-m Nordic Optical 
Telescope (NOT) of Roque de los Muchachos Observatory (La 
Palma, Spain). We employed the high-res fibre, which provides 
a resolving power of R = 67 000 in the spectral range 3600- 
7400 A. 

The SOPHIE and HARPS spectra were extracted using the 
respective pipelines. The radial velocities were then computed 
following a technique of weighted cross-correlation of the spec¬ 
tra u3ing_ajiumericalG2 sta r mask. This techn ique is described 
bv iBaranne et alJ (Il996h and lPene et al.l (l2002l) . The spectral or¬ 


Fig. 7. Radial velocity measurements and best solution for CoRoT-28. 
Offsets have been subtracted according to tbe values given in Table[6l 


ders with low signal-to-noise ratio were discarded to perform the 
computation to reduce the dispersion of the measurements. We 
eliminated the two bluest and five reddest orders for SOPHIE 
(39 in total), the ten bluest and the two reddest for HARPS (72 
in total). The cross-correlation function of CoRoT-28b shows a 
single peak with EWHM of 9.9 kms^' with SOPHIE and 7.6 
km s ' with HARPS. Also, SOPHIE and HARPS radial veloc¬ 
ities obtained with different numerical spectral masks (EO, G2 
and K5) show the same behaviour, suggesting that the radial ve¬ 
locity variations cannot be explained by a blend scenario with 
different spectral type stars. 

Eor the EIES data, we adopted the same observing strategy 
described in iBuchhave et al.l (1201 Ol) . i.e. we split each epoch ob¬ 
servation into three consecutive exposures of 1200 sec to remove 
cosmic rays hits, and we acquired long-exposed ThAr spectra in 
“sandwich-mode” to trace the RV drift of the instrument. We 
reduced the data using a customized IDE software suite, which 
includes bias subtraction, flat fielding, order tracing and extrac¬ 
tion, and wavelength calibration. We obtained radial velocity 
measurements performing a multi-order cross-correlation with 
the RV standard star HD 182572, observed with the same instru¬ 
ment set-up as the target object, and for which we adopted an 
heliocentri c radia l velocity of -100.350 km s ', as measured by 
lUdrvetalJ (119991) . 

The radial velocity measurements with the best orbital so¬ 
lution are shown in Eigure |7] assuming the period and transit 
epoch determined by the CoRoT light curve. Einally, the cross¬ 
correlation function bisector span shows no correlation with the 
radial velocities (see Eig. [8]l, reinforcing the conclusion that the 
radial velocity variations are not due to a blending effect. Radial 
velocities, their errors, and bisector span are listed in Tables IaTI 
to lA.3l 

For CoRoT-28b, we fit SOPHIE, HARPS, and PIES data to¬ 
gether while fixing the ephemeris derived from the light-curve 
analysis (see Table |6l). The rest of the orbital parameters (or¬ 
bital eccentricity, argument of periastron, radial velocity semi¬ 
amplitude, and systemic velocities) were derived using the mod¬ 
elling pr ocedure and Mar kov Chain Monte Carlo module from 
PASTIS (iDfaz et al.l2014l) . The priors used for the orbital param¬ 
eters were non-informative (uniform priors) except for the period 
and the transit epoch that were fixed. The results are summarized 
in Table |6l 

CoRoT-29b 
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Fig. 8. Bisector span as a function of radial velocities for CoRoT-28. 
We subtracted the systemic velocities in the X axis. The uncertainties 
in the bisector span is assumed to be twice the uncertainties in the ra¬ 
dial velocities. The bisector span shows no correlation with the radial 
velocities. 


Because of its faintness (V = 15.6), CoRoT-29 radial veloc¬ 
ity follow-up could only be performed with HARPS in the EGGS 
mode to improve the throughput. Compared to the HAM mode, 
which was used to follow up CoRoT-28, the EGGS mode of 
HARPS has a larger fiber (1.4" compared to 1") and no scram¬ 
bler. The spectral resolution is therefore reduced (R - 80 000) 
but the throughput is twice as large. We obtained 20 mea¬ 
surements between 19 June 2012 and 10 August 2013, show¬ 
ing that the radial velocity variations were in agreement with 
the CoRoT ephemeris and establishing the planetary nature of 
CoRoT-29b (see Pig. |9]l. The same cross-correlation technique 
with a G2 spectral type numerical mask allowed us to derive the 
radial velocities. As for CoRoT-28, the ten bluest and the two 
reddest orders were discarded. We obtained a 7.1 km s ' sin¬ 
gle peak cross-correlation function. Also, cross-correlation with 
other masks (PO and K5) show no discrepancies in the radial 
velocity variations. Moreover, the cross-correlation function bi¬ 
sector span shows no correlation with the radial velocities (see 
PigurefTOl). This suggests that the probability to have a blending 
effect is very low. The radial velocity measurements are repre¬ 
sented on Pigure 0 and are listed with their errors and bisector 
spans in Table [AA]. 

We also fit CoRoT-29b HARPS data while fixing the 
ephemeris as described above for CoRoT-28b and the results are 
also summarized in Table |6l 

3.3. Spectroscopic characterization 

To derive the stellar parameters from the spectroscopic data, 
we used the Spectroscopy Made Easy (SME) package (version 
5, Pebruary 2014). The develo pment and structure and pe rfor- 
mance of SME is described in IValenti & PiskunovI (Il996h and 
in e.g. IValenti & Pischeii (l2005h . Briefly, this code uses a grid 
of stellar models to which is fitted the observed spectrum. This 
is done by calculating a synthetized spectrum and minimizing 
the discrepancies through a non-linear, least-squares algorithm. 
The analysis is based on the latest ge neration of MARCS model 
atmospheres (Gustafsson et al.ll2008l) and the ATLAS 12 models 
(lKuruczll2005 ). The models calculate the temperature and pres¬ 
sure distribution in radiative and hydrostatic equilibrium, assum¬ 
ing a plane-parallel stellar atmosphere. 



-0.2 0 0.2 0.4 0.6 0.8 

Time (transit orbital phase) 


Fig. 9. Radial velocity measurements and best solution for CoRoT-29. 
Offsets were subtracted according to the values given in Table 



RV (km/s) 

Fig. 10. Bisector span as a function of radial velocities for CoRoT-29. 
We subtracted the systemic velocities in the X axis. The uncertainties 
in the bisector span is assumed to be twice the uncertainties in the ra¬ 
dial velocities. The bisector span shows no correlation with the radial 
velocities. 

The SME is developed so that, by matching synthetic spectra 
to observed line profiles, one can extract the information in the 
observed spectrum and search among stellar and atomic param¬ 
eters until the best fit is achieved. The SME utilizes automatic 
parameter optimization using a Levenberg-Marquardt mini¬ 
mization algorithm. Synthetic spectra are calculated on the fly, 
by a built-in spectrum synthesis code, for a set of global model 
parameters and specified spectral line data. Starting from user- 
provided initial values, synthetic spectra are computed for small 
offsets in different directions for a subset of parameters defined 
to be “free”. The model atmospheres required for this are calcu¬ 
lated through interpolating the grid of models mentioned above. 

We used as input data the high-resolution spectra taken with 
HARPS (see Section 13.21) . In the observed spectra ,we use a 
large number of spectral lines, e.g. the Balmer lines (the ex¬ 
tended wings are used to constrain Teff), Nai D lines, Mgi b, 
and Ca i (for and log g) and a large number of metal lines 
(to constrain the abundances). Por CoRoT-28, we obtained a Li i 
equivalent width of 31mA at 6 708A. 

After finding the Teftfrom the Balmer line wings (e.g. 
Ipuhrmann et alJ [l998h . we use the calibration of iBruntt et al.l 
(1201 oh to fix the micro- and macro-turbulence parameters be¬ 
fore finding other parameters. The quoted errors are 1 cr standard 
deviation errors in SME and thus a measure of how well the 
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generation of the synthetic spectrum has succeeded. We know 
(Fridlund et al., in preparation) that when we apply SME to our 
(only) known source, the Sun, we find somewhat larger errors 
such as a +100K difference in Teff, and a +0.1 dex in [Fe/H] (as 
proxy for metallicity) depending on the choice of model grid be¬ 
ing used (e.g. Atlasl2 or MARCS 2012). The mass, radius, and 
age of the star are c alculated using stellar evolution tracks from 
[Hurley et alJ (l2000h and the constraints from the analysis of the 
light curve (see SectionH]). They are reported in Table| 6 ] 

Additionally, we estimated the interstellar extinction Ay and 
spectroscopic distance d to CoRoT-28 and CoR oT-29 following 
the method described in iGandolfi et al.l (l2008l) . We simultane¬ 
ously fitted the available ExoDat, 2MASS, and WISE colours 
(Table [TJ with synthetic theoretical r nagnitudes derived by 
NextGen low-resolution model spectra dHauschildt et al.l 09991) 
having the same photospheric parameters as the stars. We ex¬ 
cluded the Wj, and W 4 WISE magnitudes because of the low S/N 
ratio. Assuming a normal extinction {Ry = 3.1), we found that 
CoRoT-28 and CoRoT-29 are at 560 ± 70 pc and 765 ± 50 pc 
from the Sun, and suffer an interstellar extinction of Ay - 
0.75 ± 0.20mag and Ay = 0.85 ± 0.15mag, respectively. The 
dereddened B-V excess of the stars is 0.98 + 0.06 and 0.87 ± 0.03 
mag, respectively. 

3.4. Spectral typing 

A spectral type was derived for CoRoT-28 from a spectrum taken 
on 24 July 2012, with the low-resolution spectrograph (^ » 
1000) at the Nasmyth focus of the 2m telescope in Tautenburg. 
A slit width of l"and the V200 grism were used covering the 
wavelength range 360-935 nm. The signal-to-noise of the spec¬ 
trum is about 180. The sp ectrum was reduced and extracted with 
IRAF (iTodvIl 1986111^ in the standard way, including subtrac¬ 
tion of the bias and the background as well as flat-held and ex¬ 
tinction correction. Wavelengths were calibrated using sky lines 
in the long-slit spectra and exposures of He and Rr gas discharge 
lamps. The spectral type of CoRoT-28 was obt ained following 
IGandolfi et al.l (l2008l) and lSebastian et all d2012h . Briefly, we fit¬ 
ted the observed lo w-resolution spectru m using template spectra 
from the library bv IValdes et al.l ( 20041) in combinati on with the 
photospheric parameters derived bv I Wu et alJ (1201 ih . Although 
a best match was attained for G8/G9V templates, we also found 
matching templates with spectral classes G5-K0 and luminos¬ 
ity classes IV and III. The luminosity class is not well defined 
by the low-resolution spectrum. For the reasons explained above 
and considering the transit parameters, we consider that the star 
is slightly evolved and classify it as G8/9IV. For CoRoT-29, 
considering the effective temperature , we ch oose a spectral type 
KOV following iKenvon & H^tmamild 19951) . 

4. Planetary parameters 

4.1. CoRoT-28b 

For the transit analysis and final planetary parameter determi¬ 
nation of CoRoT-28b, we used the light curve that was filtered 
for stellar variability, as explained in Section |2]2] and the infor¬ 
mation from the spectral analysis described in Section [3 Since 
no significant transit timing variations (TTVs) were found, we 
folded the light curve using the ephemeris in Table| 6 ] The folded 
light curve was then mo delled with the Transit Light Curve 
Model code, described in ICsizmadia et al.l d201 ll) . Each chro¬ 
matic light curve (red, green, blue, and white, see Section HD 
was treated in the same way. For the transit light-curve fit, we 


used the iMandel & Agoll d2002h model. The optimization pro¬ 
cess consists of two steps. First, the best solu tion is found with a 
Harmony Search analysis genetic algorithm dGeem et al.ll2001ll . 
Second, the uncertainties of the parameters are obtained via sim¬ 
ulated annealing around the best solution found by the genetic 
algorithm. 

The free parameters are the scaled semi-major axis {ajR^, 
where a is the semi-major axis and R^ is the stellar radius), the 
planet-to-stellar radius ratio (k = Rp/Rs, Rp is the planetary ra¬ 
dius), the impact parameter b, where b - a cos i Vl - e^/Rs(l + 
e sin a>), i is the inclination, e is the eccentricity, and w is the ar¬ 
gument of the periastron. Values of e and a> are known from the 
radial velocity analysis, and their uncertainties were propagated 
allowing them to vary between their ±lcr uncertainties during 
the optimization procedure. The epoch could also vary within 
the un certainties reported in Table |6] Following ICsizmadia et al.l 
(l2013l) . limb darkening c oefficients were fitte d as free parame¬ 
ters. W e took the results of iBrown et al.l (l2001l) and iPal & KocsisI 
(l2008h into account, fitting the combinations «+ = «„ + Uj, and 
u^ - Ua - Ub rather than the linear and quadratic coefficients of 
Ua and Ub individually. A quadratic limb darkening law provided 
a reasonable fit to the data. 

The fit is shown in Figure [TD and the results can be seen 
in Table H The fits obtained in different colours are in agree¬ 
ment with each other. Note that because of the non-axisymmetric 
chromatic PSF of CoRoT, the contamination factors are not ad¬ 
ditive, so contamination in white is not equal to the sum of con¬ 
taminations in different colours. The white light curve results 
were used to establish the planetary parameters given in Ta¬ 
ble |6] because they have the smallest uncertainties because of 
the higher SNR. 

4.2. CoRoT-29b 

As described in Section l231 the light curve of CoRoT-29b has a 
non-flat, mid-transit shape. 

The CoRoT data sets sampled at 512s and 32s cadence were 
treated separately because of the background correction issue de¬ 
scribed in SectionHD The slope of the bottom part of the transit 
is non-zero with 95% confidence level in the 32s sampled data 
by CoRoT, and only non-zero with Icr significance in the 512s 
sampled data set. It is also non-zero with 99% confidence level in 
the ground-based data from FTN, and non-zero with 95% confi¬ 
dence level in the IAC80 data (see Section lXTT l. The slopes in the 
CoRoT data and in the ground-based data are compatible with 
each other within Icr uncertainties, though the data from ground 
and from space were taken in slightly different wavelengths and 
suffer from different kind of systematic effects. 

We have not been able to identify any systematic resid¬ 
ual in the CoRoT data, which could account for the asym¬ 
metry of the transit; neither identified any other target in the 
CoRoT data set with a similar asymmetry. Ground-based ob¬ 
servations are affected by correlated noise, which c an produce 
asymmetries in the transit shape (see, for exampl e, iHebb et ^ 
I 2 OIOI : iDittmann et al.l 120101: iMancini et al.l l2013l) . But in this 
case, despite the inherent limitations of the ground-based pho¬ 
tometric observations that we present (see Section lTTI) . we con¬ 
sider that the asymmetry observed by CoRoT is confirmed by 
the ground-based observations, and therefore it could have an 
astrophysical origin. 

Possible interpretations for the asymmetry of the transit are 
given in the next section (14.31) . if the 2-sigma detection is as¬ 
sumed to be significant. The model that better reproduces the 
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Table 2. Parameters of the fit to the transit light curve of CoRoT-28b. 


Parameter 

in blue 

Value and uncertainty 
in green in red 

in white 

aIRs 

7.1 + 1.7 

6.7+ 2.7 

7.40 + 0.68 

7.29 + 0.16 

k = Rp/Rs 

0.0614 + 0.0073 

0.0600 + 0.0107 

0.0581 +0.0033 

0.0551+0.0004 

b 

0.43 ± 0.36 

0.54 + 0.47 

0.30 + 0.26 

0.24 + 0.09 

u+ 

0.42 ± 0.32 

0.71+0.40 

0.43 + 0.39 

0.78 + 0.11 


0.32 ± 0.62 

0.18 + 0.59 

0.40 + 0.58 

0.37 + 0.22 

contamination [%] 

6.2+ 3.0 

6.4+ 3.0 

2+1 

3 + 1 

^red 

1.24 

1.28 

1.20 

1.32 




0.00 

Orbital phase 


white 


0) 

cc 
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Fig. 11. The folded light curve of CoRoT-28b together with the best fit in the different chromatic light curves observed by CoRoT. Blue points are 
the observed data; the orange curves represent the individual fits. 


data is that including the gravity darkening produced by an 
oblate star (see Section |4. 3. 11 1. We have taken those values as 
the reference planetary parameters (see Table|6]l. 

4.3. Interpretations of the transit light curve of CoRoT-29b 
4.3.1. Stellar oblateness 

We could interpret the shape of the light curve as the result of 
the transit of the planet above a non-spherical stellar surface. A 
spherical surface becomes ellipsoidal, in a first approximation, 
when it is deformed by its rotation, by its internal mass distribu¬ 
tion, or by the tidal influence of a massive companion. Conse¬ 
quently, the stellar poles are closer to the centre of the star than 
the equator, thus the absorption rate of photons from the stellar 
interior is different. Additionally, the effective gravitational ac¬ 
celeration is higher at the poles than at the equator, where the 


the centrifugal force reaches its maximum. The decrease of the 
effective gravitational acceleration at the equator affects the stel¬ 
lar atmospheric scale height, which in turn changes the stellar 
temperature layering. All these effects cause a pole-to-equatorial 
variation in the stellar flux, which is called gravity darkening, 
because this temperature latitudinal variation is characterized by 
the effective surface gravity g at a certain latitude b and with a 
gravity darkening exponent fi. 


T(b) = Tpoie - 

\^pole 

The gravity d arkenin ! 


( 1 ) 


effec t was studied by von Zeinell (Il924t) : 


iLucvl (1 19671) : IClared ll999l) . and IClaretl (120121) from theoreti- 
cal point of view. In the context of tra nsiting e xtraso lar plan¬ 
ets, this effect was studied in depth by iBarnesI (l2009l) . Impor¬ 
tant observational results, in eclipsing binary s tars and transit- 
ing extrasolar planets, have been published by iRafert & Twigd 
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( 1980l): Nakamura & Kitamura ( 19921): Diurasevic et al.l (2003): 
Szabo et al.| (|20n'l)r Morris et al.l (1201 3h : Barnes et alT ( 2013 ): 
Zhou & Huand ( 2013l) . 

As is shown hereafter, we successfully explained the ob¬ 
served transit light-curve shape of CoRoT-29b with the assump¬ 
tion that the asymmetry of the light curve is caused by gravity 
darkening effect. We modelled this effect with our own code, 
which we briefly describe next, and which will be describe in 
detail in a future publication (Csizmadia, in prep.). The opti¬ 
mization method is identical to the one described before (see 
Section f4.1b . 

For this ht, we used the spectroscopic constraints for the 
vsin/ and TefTfrom the spectroscopic analysis of the star. We 
allowed these parameters to vary between 3.5 + 0.5 kms“* and 
5260 + 70K, respectively. Here Teff is the surface temperature of 
the star averaged for the visible hemisphere. 

The stellar shape is dehned by V ^constant equipotential sur- 
faces. We use the quadrupolar approach (e.g. IZahn et ^ 


instead of the simple ellipsoidal approximation of BarneJ i 20091) . 


20101 ) 


V = 


GM, 


R(b) 


I-J 2 


^s.eqV 

ml 


Flisin b) 


-D.m?R^{b)co^ b. (2) 


This equation can be rewritten in a dimensionless form as 


T' = 


VR 


s,eq 


GM, 


r(b) \ r(b)^ 


"2 (fl/Rs.eq)^ 


(3) 


Here we used the following notations: 

- Rs,eq, equatorial radius of the star, 

- G, gravitational constant, 

- Ms, mass of the star, 

- r(b) - R{b)/Rs^eq, the normalized radius of the star at latitude 

b, 

- w - Poih/Pmt, a rotational parameter, namely the ratio of the 
orbital period of the planet and the rotational period of the 
star, 

- a/Rs.eq, the scaled semi-major axis where a is the semi-major 
axis of the planet, 

- J 2 , the second (or quadrupole’s) gravitational momentum of 
the star, 

- P 2 , the second Legendre-polynomial. 

The free parameters are J 2 ', the gravity darkening exponent 
yS; the planet to stellar radius ratio k - Rp/Rs', the s tellar inclina- 
tion angle 4 (the obliquity of the star, as dehned in lBarnej|2009l 
is i/j = 90-4); , t he sky projected angle of the stellar spin axis 

(T in the notation of iFabrvckv & Winnll2009l) : ip , the inclination 
of the planetary orbit; the limb darkening coefficients m+ and m_; 
the epoch of the transit Tq; the temperature of the star at the pole 
Tpoiei and the scaled semi-major axis ci/Rtj gq(see also Figs. IBT] 
and [T2]) . The rotational parameter w was calculated iteratively 
from the known orbital period and from v sin /, assuming a cer¬ 
tain average radius of the star. We chose a stellar radius of 1 
solar radius as initial value. After the hrst ht, we calculated the 

density parameter M'^^/Rs eq = |^3;7r(a/Rs,eq) /GForb^] for the 
star. We searched stellar models by iHurlev et al.l (l2000l) . which 
provide this density parameter at the effective surface tempera¬ 
ture of the star (within the uncertainties), and this yielded a new 
estimate for the stellar radius. With this new stellar radius, we 
repeated the ht until convergence (see Fig. ED. 


Table 3. Parameters of the ht to the transit light curve of CoRoT-29b 
with a model accounting for gravity darkening. 


Fitted parameters 

Value and uncertainty 

H 

0.028 + 0.019 

Gravity darkening exponent [i 

0.89 + 0.25 

k = Rp/Rs 

0.1028 + 0.0043 

4 [°]“ 

46+ 19 

a [°]'’ 

256 + 40 

4 [°] 

87.3 + 2.7 

u+ 

0.62 + 0.15 

U- 

0.58 + 0.14 

7’o[s] 

-200+31 

Tpole [K] 

5 341 + 266 

^/^s,eq 

9.22 + 0.19 

contamination factor [%] 

54 + 4 

Derived parameters 


f ~ ^s,pole/^s,eq 

0.94 + 0.02 

w = Foi-bZ-Prot 

0.26 + 0.06 


71 

2 c 

/^red 

1.4 

BIC 

120 


Notes. The obliquity of the star, as dehned in iBarne i ll2009f) . is ip = 
90 - 4 = 47°. The sky projected ang le of the stellar spin axis; A in 
the notation of iFabrvckv & Winnl (l2009h . Sixty-two htted points and 
12 free parameters. 


The stellar parameters dehne the value of T'. A Newton- 
Raphson iterative process yielded the value of r and the actual 
value of the surface gravity g = W for every latitude and lon¬ 
gitude. Then we calculated the effective surface temperature as 
well as the intensity assuming a simple black-body radiation for 
the star. A numerical integration procedure provided the unob¬ 
scured Hux of the star with a freely chosen gravity darkening 
exponent and limb darkening coefficients. As is well known, the 
limb darkening coefficients are functions of the effective stellar 
surface temperature. If the temperature varies over the surface, 
different limb darkening coefficients will be valid at each sur¬ 
face point. Since the theoretical limb darkening tables are not 
always in agreement with each other and they have not been val¬ 
idated observationally, it is hard to accept that an interpolation 
of these tables will provide the appropriate description of the 
limb darkening coefficients for every surface point. In addition, 
a rotating stellar atmosphere can be different from a static at¬ 
mosphere, whereas most limb darkening tables use non-rotating 
stellar atmospheric models. Therefore, we left the limb darken¬ 
ing coefficients free and used a quadratic law hoping that this 
will be a better approximation of reality in this particular case. 
For furthe r details about th e diffic ulties of limb darkening han¬ 
dling, see ICsizmadia et al.l (l2013l) . Our code also includes the 
rotational beaming effect. 

The planet’s position was projected on the sky and when 
it crossed the apparent stellar disc, we calculated the amount 
of the blocked stellar e mission in a small-planet approximation 
(iMandel & Agolll2002l) . Then we adjusted the free parameters 
until we had a convergence in the ;^f^-values, using the optimiza¬ 
tion methods described before ('Section [4.1b . The results of our 
modelling are shown in Table |3] and in graphical form in Fig¬ 
ure [T3] The hnal planetary parameters are also shown in Table |6] 
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Fig. 12. Angle convention used in the gravity darkening model. 
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Fig. 13. The binned and folded light curve of CoRoT-29b together with 
the best fit accounting for gravity darkening fit. See text. 


Table 4. Parameters of the fit to the transit light curve of CoRoT-29b us¬ 
ing a spherically symmetric model for the star. We discarded this model 
as a correct explanation of the measured light curve, as explained in the 
text. 


Parameter 

Value and uncertainty 

aIRs 

10.49 + 0.72 

Rp/Rs 

0.0838 + 0.0070 

b 

0.08 + 0.27 

u+ 

0.89 + 0.33 

U- 

-0.49 ± 0.38 


103 

2 a 
^red 

1.87 

BIC 

132 


Notes. Sixty-two fitted points and seven free parameters. 



Fig. 14. The binned and folded light curve of CoRoT-29b together with 
the best fit using a spherically symmetric model of the star. We dis¬ 
carded this model as a correct explanation of the measured light curve 
as explained in the text. 


4.3.2. Spherical star model 

If we decide to ignore the non-flat shape of the transit and try to 
model the light curve of CoRoT-29b with the same techniques 
used in the case of CoRoT-28b (see Section 14. Il l, we obtain a 
fit solution significantly worse than using the more sophisticated 
analysis including gravitational darkening. The results are sum¬ 
marized in Table ID and in graphical form in Fig. [14] Notice the 
remarkable distribution of the binned points relative to the fit. 
After the ingress, the observed points are systematically below 
the fit, and later, during the mid-eclipse and just before the egress 
phase they are systematically over the best match of the spheri¬ 
cal star fit. We conclude that we cannot ignore the asymmetry of 
the light curve in the observed data set. 


4.3.3. Spot-occulting model 

Stellar spots occulted by transiting planets can have a significant 
impact in the measured transit light curve. This phenomenon has 
been widely studied i n the literature and has been con firmed by 
several observations ( Silva! 1200^ Wolt er et al.ll2009l and refer¬ 
ences therein, ISanchi^Oieda et al. 1201 2i etcA We require from 
a spot-occulting solution to reproduce the observed light curve 
with a stellar spot coverage, which is consistent with our knowl¬ 
edge of the stellar properties. Typically, the parameters used for 


this models are the stellar temperature and rotational period and 
the spot temperature, size, lifetime, and distribution. 

In particular, the typical lifetime of a spot is ~ 2 weeks for 
the Sun, but other stars have both s horter and longer li fetimes, 
up to several months in some cases dOueloz et alJl200lh . In our 
case, comparing the identical CoRoT and the FTN observations 
(see Section 13.11 and Fig. O, we require a spot lifetime larger 
than one year. 

Secondly, as all planetary transits observed by CoRoT in 
2011 are indistinguishable with the current precision level, and 
they are also indistinguishable from FTN data from 2012 and 
IAC80 ground-based observations from 2014, we must conclude 
that the planet transits every 2.85 days above the same spot. This 
requires an exact synchronization between the stellar rotational 
period and the planetary orbital period, with exactly the same 
phase in a baseline larger than three years. 

On the one hand, if the observed asymmetry is related to an 
equatorial or mid-latitude stellar spot, then the stellar rotational 
period should be equal to the orbital period of the planet, other¬ 
wise we cannot observe the spot always in phase with the plane¬ 
tary transits. However, we do not see a signature of such a spot as 
a modulation of the light curve in phase with the stellar rotational 
period. Moreover, we require the star to have a rotational period 
identical to the orbital period, which is not consistent with the 
stellar parameters derived from spectroscopy (in particular, with 
the vsin/)- 
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On the other hand, if the spot is polar and the rotational axis 
of the star is perpendicular to the orbital plane of the planet, then 
we could not observe an asymmetric light curve. A spot-caused 
asymmetry requires an off-position of the spot from the apparent 
stellar disc. This scenario has the advantage that the signature 
of the spot is diluted in the light-curve signal, and we could not 
infer from the modulation of the light curve the stellar rotational 
period. Additionally, if the spot is centred on the pole and the 
stellar rotation axis has the right orientation, the planet could al¬ 
ways cross the same location of the spotted region. This requires 
no prior assumptions on the stellar rotation period (see Fig. [BJ. 

If we accept the polar spot scenario, however, which a priori 
cannot be ruled out as easily as the equatorial spot model, then 
we must require that the polar spot does not change its apparent 
position during the baseline of the observations. In this model, 
the rotational axis of the star is very different from the angu¬ 
lar momentum of the vector of the planet, or in other words, 
the planetary orbit is misaligned with the equatorial plane of the 
star. The spot cannot be at the centre of the apparent stellar disc 
because of the asymmetry of the light curve. According to the 
light curve, the planet crosses the unspotted surface of the star 
hrst, producing a deeper transit, and later the spotted surface, 
when the lower brightness temperature of the spot produces a 
shallower transit. The apparent size of the spot should be much 
larger than the apparent size of the planet because we observe a 
continuous change in the transit depth and not a step-like event. 

Polar spots in active stars are a common phenomenon and 
in some cases they can survive for long intervals of time. Or at 
least, active regions continuously producing spots at similar lati¬ 
tudes have been obse rved for sever a l years in some y oung active 
stars like V410 Tau (lHatzesiri995t iRice et al.ll201 Ih . But these 
polar spots are on very active, young stars that are rapidly rotat¬ 
ing (V410Tau has a vsin/of 75kms“'). Long-lived polar spots 
have never been seen, to our knowledge, on more slowly rotating 
stars like CoRoT-29. Moreover, these spots should produce sig¬ 
natures on the stellar spectrum (Ca ii emission, distortion in the 
line profiles, flat-bottomed line profiles), which are not observed 
in our data set (see Section[3]). 

Although it is unlikely that a spot can survive for such a long 
period without evolving on the surface of a main-sequence star, 
we used this assumption in our modelling with a spherical star, 
a spherical planet, and a circular spot on the surface of the star, 
with an arbitrarily oriented stellar rotational axis and arbitrary 
initial spot properties. We use a quadratic limb darkening law, as 
before. The results are given in Table |5] and shown in graphical 
form in Fig. [161 

We obtain a polar spot with a tilted stellar rotational axis. 
This tilt means here that the projected rotational axis of the star 
and the projected angular momentum of the planet do not coin¬ 
cide. The planet is on ’pole-on’ orbit because it orbits the star 
in such a way that its projected orbit crosses the polar region of 
the star. The mutual inclination between the stellar equator and 
planetary orbital plane is practically 90 degrees. Note that such 
an orbit is very stable against perturbations. 

From the measured v sin i and from the fitted stellar inclina¬ 
tion for the spotted case (assuming ~ 0.9Ro), we find that 
the stellar rotational period should be between 0.7-8.4 days. The 
spectroscopically measured vsini ~ 3.5 kms“' would yield 
Veq ~ 10 kms“^ with 4 » 20°, where the inclination of stel¬ 
lar rotational axis comes from the spot fit. 

The comparison of the Bayesian Information Criterion (BIC) 
of the spot fit and the gravity darkening fit makes a clear choice 
which fit is better. The spot fit yielded BIC - 133, whereas the 
gravity darkening fit produces BIC - 120. This BIC difference 


Table 5. Parameters of the fit to the transit light curve of CoRoT-29b 
using a spherically symmetric model for the star and a stellar spot. The 
parameters 4pot and dspo, are the astrographic latitude, longitude, 
and diameter of the spot in degrees. The stellar angles are chosen in the 
same way as in Fig. [HI We discarded this model as a correct explanation 
of the measured light curve, as explained in the text. 


Parameter 

Value and uncertainty 

aIRs 

10.32 + 0.24 

Rp/Rs 

0.0948 + 0.0031 

ip 

90.7 + 0.8 

u+ 

0.90 + 0.19 

U- 

-0.26 + 0.28 

T, [K] (fixed) 

5260 

Tspot [K] 

1600 + 400 

^spot [°] 

-56 + 21 

kpol [°] 

209 + 45 

r^spot [°] 

39+13 

4 [°] 

-42 + 20 

a n 

211+39 

contamination factor [%] 

45 + 3 


84 

2 a 
/Vred 

1.75 

BIC 

133 


Notes. Sixty fitted points and 12 free parameters. 


corresponds to a Bayes factor of around 665 for the gravity dark¬ 
ening model (i.e. if the prior probabilities for each model are the 
same, then the gravity darkening is 665 times more probable than 
the spot model). Moreover, to explain the observational evidence 
we are obliged to put requirements on the size and temporal evo¬ 
lution of the spot, which are completely ad hoc. We consider that 
there is enough evidence to discard the spotted star as a reason¬ 
able interpretation of the data. 

An alternative to the single polar spot would be a band of 
spots with short lifetimes, but continuously appearing during the 
baseline of our observations in the same latitudes, relaxing one 
of the requirements of the scenario. One possibility is that if the 
spots were conveniently appearing at different longitudes, the 
condition of the synchronous rotation could also be relaxed. In¬ 
cluding more spots will improve the modelling results at the cost 
of increasing the number of free parameters. Although plausi¬ 
ble, it is finally also an ad hoc solution for the problem. This 
hypothesis could be falsified, however, if observations of the 
Rossiter-McLaughlin effect would show a small spin-orbit an¬ 
gle. 


4.3.4. Planetary oblateness 

A possible wind-driven oblated shaped of the planet can cause 
asymmetries in t he light curve, but t he observable am plitude is 
about 100 ppm (iBarnes et al.ll200^ IZhu et al1l2014ll . which is 
about 200 times smaller than the observed effect in our case, 
therefore we reject this explanation. 

Other examples of star-planet interaction do not provide 
satisfactory answers in our case because they have typically 
different amplitudes and timescales (i.e. iJackson et al.l 1201^ 

lEsteves et al.ll2013tlFaigler & Mazehll2015h. 
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Fig. 15. Visualization of the CoRoT-29 system according to the spot 
model. The scale is in stellar radii. The orange area is the spotted area, 
while the horizontal line shows the sky-projected planetary orbit. The 
positions of the planet at different orbital phases are shown; the size of 
the circles correspond to the planetary radius in this scale. The pole of 
the star is denoted by a small green cross. 



Fig. 16. The binned and folded light curve of CoRoT-29b together with 
the red line being the best fit using a spherically symmetric model for 
the star and a stellar spot. For comparison, the solid green curve is be¬ 
cause of the unspotted light curve model (we set rspot= T^in the fitting 
routine). We discarded this model as a correct explanation of the mea¬ 
sured light curve, as explained in the text. 


4.3.5. Other scenarios 

A disc around the star could produce an asymmetric transit, but 
there are no signs of such disc in the light curve or in the infrared 
emission of the star. A moon around the planet could produce an 
asymmetric transit, but dynamically there is no room for a moon 
around a planet so close to its star. A second planet in the system 
transiting simultaneously could produce an asymmetric transit, 
but not every time that the transit is observed, and there are no 
signs of this planet in the radial velocity. We conclude that none 
of these scenarios provides a satisfactory explanation of the data. 


5. Discussion 

5. 1. Stellar properties 

CoRoT-28 

CoRoT-28 is a G8/9IV evolved star with a radius of ap¬ 
proximately 1.8 solar radii. With an age of approximately 12 


Gyr, we expect from gyrochronology a rotational period around 
70 or 80 days and a vsin/around 1.2kms“* or smaller (iBarnesI 
I2OO7I) . The V sin i value from spectroscopy is almost three times 
larger, giving an expected rotational period of about 30 days. 
The autocorrelation f unction (ACF) of the light curve (see 
iMcOuillan et al.ll2013h does not show any evident periodicity at 
all, but if there is any sign, it is for periods larger than the 85 
days observed by CoRoT{see Fig. [nil. 

There is an inconsistency between the v sin / determined by 
spectroscopy, which is consistent with a relatively fast rota¬ 
tional period, and the stellar parameters derived from the mod¬ 
els, which suggest an old star with a slow rotation rate. From 
our analysis of the light curve we cannot decide between the 
two scenarios. The v sin/has a relatively large uncertainty, but 
seems to be consistent with the higher rotation rate. It is pos¬ 
sible that the stellar models fail to reproduce the actual age of 
the star. The range proposed by the models is 12.0 +1.5 Gyr, 
close to the old end of the distribution of stellar populations, 
which is consistent with the evolved stellar radius from the mod¬ 
els (1.78 + 0.11 Rjfor an effective temperature of 5 150 + 100 
K). This range, however, is perhaps not as consistent with the 
metallicity [Fe/H] = 0.15 + 0.05 dex, which is quite a robust 
observational result. 

There is also an apparent inconsistency between the age and 
the Li I abundance in the at mosphere of the s t ar. We estimated 
logA(L/) using Table 2 of ISoderblom et aP (Il993h given the 
equivalent width of 31mA. For a dwarf with an effective temper¬ 
ature of 5 150K, we obtained logN(Li) - 1.4. This is in line with 
an age of a few hundred Myrs according to ISestito & Randichl 
(l2005h . However, stars of this effective temperature with ages 
of few hundred Myrs have already reached the main sequence, 
which is not compatible with the surface gravity measured from 
spectroscopy or with the transit duration, which indicates that 
the star has a radius larger than the Sun. Moreover, we would 
expect a significant degree of stellar variability and a faster rota¬ 
tion period, which is also not supported by the observations. The 
abundance of lithium and its relationship with age, especially 
for stars hosting planets (see, for example, t he recent discus¬ 
sions by iLiu et al.ll2014l: lAdamow et al.ll2()l^ . is still debated. 
Therefore we do not consider lithium as a reliable age indica¬ 
tor for CoRoT-28. On the other hand, it is known that a small 
amount of Li i can also be present in evolved stars, for example 
lithium can be regenerated by the Cameron-Fowler mechanism 
(ICameron & Fowle3ll971h . 

We consider that observational evidence speaks in the case 
of CoRoT-28 in favour of an evolved star with a higher lithium 
content and v sin / value than what would be expected for its age 
and evolutionary status. 

CoRoT-29 

CoRoT-29 is a KO main-sequence star with a radius of ap¬ 
proximately 0.9 solar radii. Its age is not well constrained from 
the analysis of the stellar models and the autocorrelation func¬ 
tion of the light curve (see Fig. [TST i does not help to constrain 
the rotational period. In any case, the analysis of the ACF in the 
case of CoRoT-29 can be debated, as ~ 50% of the flux within 
the mask comes from a contaminant. 

It is however interesting that the analysis of the light curve 
(Section I4.3l l shows that the orbital plane of the planet could 
be misaligned with respect to the equatorial plane of the star. 
The gravity darkening model allows us to measure the relative 
inclination of the spin axis of the star, a parameter that typically 
cannot be constrained with photometric measurements. 
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Table 6. Planetary and stellar parameters. 



CoRoT-28b 

CoRoT-29b 

Ephemeris 



Planet orbital period P [days] 

5.208 51 +0.000 38 

2.850570 + 0.000006 

Transit epoch [HJD-2450000] 

Transit duration [h] 

Results from radial velocity observations 

5 755.846 6+0.002 9 
5.70 + 0.09 

5753.115+0.001 

2.19 + 0.04 

Orbital eccentricity e 

^■^^/_0O38 

0.082 + 0.081 

Argument of periastron u [°] 

126+'f 


Radial velocity semi-amplitude K [ms'] 

56.4 + 4.9 

125 + 17 

Systemic velocity Vr (HARPS) [kms^*] 

76.7081+0.0051 

64.1907+0.0012 

O-C residuals (HARPS) [ms '] 

^ - 4.8 

11 + 11 

Systemic velocity Vr (Sophie) [kms^*] 

O-C residuals (Sophie) [ms“*] 

76.7549+0.0082 
32.8 + 6.4 


Systemic velocity Vr (FIES) [kms^*] 

76.8642+0.0074 


O-C residuals (FIES) [ms“'] 

Fitted transit parameters 

8.5+ 8.3 


Scaled semi-major axis a/Rs,eq 

Radius ratio k = Rp/R^ 

7.29 + 0.16 
0.0551+0.0004 

9.22 + 0.19 
0.1028+0.0043 

Quadratic limb darkening coefficients® u+ 

0.78 + 0.11 

0.62 + 0.15 


0.37+0.22 

0.58+0.14 

Impact parameter'’ b 

Orbital inclination ip [deg] 

0.24 + 0.09 

87.3 + 2.7 

Stellar rotational axis inclination f [°] 

- 

46+19 

Stellar rotational axis orientation O, [°] 

- 

256 + 40 

Gravity darkening exponent f 

J 2 for the star 

— 

0.89 + 0.25 
0.028+0.019 

Contamination factor [%] 

Deduced transit parameters 

3 + 1 

54 + 4 

Stellar radius ratio / = Rs,poie/^s,eq 

1 (fixed) 

0.94 + 0.02 

Mf^^IRs [solar units]® 

0.58+0.02 

1.09 + 0.02 

Stellar density ps [kgm“^] 

270 + 6 

1 825 + 38 

Impact parameter'’ b 

Orbital inclination ip [deg] 

Spectroscopic parameters 

88.1+0.8 

0.43 + 0.44 

Effective temperature Teff [K] 

Stellar surface gravity log g [cgs] 

Metallicity [Fe/H] [dex] 

5 150+100 

3.6+ 0.2 

0.15 + 0.10 

5 260+100 

4.3+0.2 

0.20 + 0.10 

Stellar rotational velocity v sin i [ km s'] 

3.0+ 0.5 

3.5+0.5 

Microturbulence velocity Fmic [ km s“'] 

0.9 + 0.1 

0.9 + 0.1 

Macroturbulence velocity Vmac [kms^'] 

3.3+0.3 

1.2 + 0.3 

Lithium equivalent width [mA] 31 

Spectral type G8/9IV 

Stellar and planetary physical parameters from combined analysis 

no Li at 6708A 

KO V 

Star mass [Mq]® 

1.01 +0.14 

0.97 + 0.14 

Star radius [Rq]® 

1.78+0.11 

0.90 + 0.12 

Stellar surface gravity log g [cgs] 

Age of the star t [Gyr] 

3.94 + 0.12 
12.0+1.5 

4.52 + 0.19 

1-8 

Distance of the system [pc] 

560 + 70 

765 + 50 

Interstellar extinction Ay [mag] 

0.75+0.20 

0.85 + 0.15 

Orbital semi-major axis (Kepler’s law) a [AU] 

0.0590 + 0.0027 

0.0390 + 0.0020 

Orbital semi-major axis (light curve) a [AU] 

0.0603+0.0050 

0.0386 + 0.0059 

Planet mass Mp [Mjup]®' 

0.484 + 0.087 

0.85 + 0.20 

Planet radius Rp [Rjup]®' 

0.955+0.066 

0.90 + 0.16 

Planet density pp [ g cm“^] 

Planet surface gravity logg [cgs] 

0.69 + 0.27 

3.12 + 0.14 

1.45 + 0.74 

3.42 + 0.19 


Notes. /(/j)//(l) = 1 - /r + UiiH + ui,( \ -yu)^, where /(I) is the specific intensity at the centre of the disc and ji = cosy, y being the angle between 
the surface normal and the line of sight; «+ = u „ + Uh and = Wa - nj. ® h = a cos i Vl - + e sinoi) Radius and mass of Sun taken as 

695 500 krn and 1.9891 x 10'^° kg, respectively (lLanglll91^ . Radius and mass of Jupiter taken as 71 492 km and 1.8992 x 10^’ kg, respectively 
( lLanalll999h . 
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Fig. 17. Autocorrelation function of the light curve of CoRoT-28. 

The alignment between the vector normal to the planetary or¬ 
bital plane and the stellar spin axis is considered one of the key 
observables constraini ng planetary for mation and evolution (for 
a recent discussion, see lDawsonl2014l and references therein). It 
has been proposed that the alignment, or misalignment, is related 
to the dynamical interaction between the star and the planet, 
which in turn depends on the internal distribution of matter in¬ 
side the star. Cold, convective stars hav e different dissipatio n 
timescales than hot, radiative stars (see lAlbrecht et al.l 12012h. 
CoRoT-29 with a Tefr of 5 260K and with a misalignment angl^ 
of Qs = 301 + 30°would be right between the out liers of the dis- 
tribution HAT-P-11 and HD 80606 in Figure 20 of lAlbrecht et al.l 
( I 2 OI 2 I) . It is however worth mentioning that the tidal dissipation 
timescale of CoRoT-29b is 1000 times faster than for HAT-P-11, 
using equation 2 from the same paper. As mentioned, the dissi¬ 
pation timescale strongly depends on the internal structure of the 
star, which deserves further attention in the case of CoRoT-29. 

The shape of the transit light curve of CoRoT-29b has a 
strong contribution from the gravity darkening, or ders of mag¬ 
nitud e larger than in othe r examples found (i.e. ISzabo et^ 
I 2 OI Ul^ou & Huang||2013h . The contributions are only compa- 
rable in the case o f the pre-main-sequence star PTFO 8-8695 
dBarnes et al.ll2013l) . However, the origin of the gravity darken¬ 
ing effect in PTFO 8-8695 is the fast rotation of the star, a M 
dwarf with a radius of 1.4Ro, a rotational period of 10 hours, and 
a Vsin/of 80kms *. For CoRoT-29b, the origin of the extreme 
pole to equatorial radius difference is not as easily understood 
(see Fig. [13. The spectroscopically measured v sin / and the re¬ 
sulted stellar inclination (43°) yield Veq = 5.1 + 1.8 kms^' and 
that corresponds to a stellar rotational period of 4-13 days. 

To lowest order, and in the assumption of uniform rotation, 
the stellar shape is defined by the surface of constant total po¬ 
tential as defined in Eq.|2] It shows the balance between the con¬ 
tribution to the local gravitational potential of the stellar rotation 
Drot and the quadrupole moment term dominated by 72- Actu¬ 
ally, we should have included in the balance the contribution 
from the tidal distortion caused by the planet, and account for 
the misalignment between the rotation axis of the star and the 
angular momentum of the orbit. But when comparing the tidal 
and rotational contributi ons to the distortion (see equa tion 3 of 
iRagozzine & Wolil2009l or A.21 of iLeconte et al.ll201 Ih . it turns 
out to be a couple orders of magnitude less important, and so we 


* fls is the sky projected an gle of the stellar s pin axis, A in the notat ion 
of iFabrvckv & Wim i2009l) . In the notation of lAlbrecht et all i2012l) . it 
would be /I = 180 - Q.^. 


neglect it. Note also that in equation [2] the contribution from ro¬ 
tation is negligible compared to the contribution of the J 2 term. 
This means that in our modelling the rotation rate is constrained 
mainly by the v sin i value. The model can measure the deforma¬ 
tion and the orientation of the star, but not its rotation rate. The 
results of Section 14.21 show that for CoRoT-29, the quadrupole 
moment has a contribution 2 orders of magnitude larger than ro¬ 
tation. Actually, the quadrupole moment of CoRoT-29 measured 
from the gravity darkening is far too large, J 2 - 0.028 ± 0.019, 
compared to the Sun, 7® = (1-7 + 0.4) ■ 10“^ ( Lang|[l99^ . Other 
authors have claimed large Jj values for fast-rotating stars, such 
as WASP-33 {J 2 - 3.8 ■ 10 '^. Ilori^l201 111 , but not as large as for 
CoRoT-29. 

The value of the quadrupolar moment depends on the dis¬ 
tribution of matter inside the star. When the star is in uni¬ 
form rotation, this quadrupolar moment is closely related to the 
Love numbeiQ ^2 in the linear approximation. Indeed, once ^2 
is known, the external potential and the shape that a body ex¬ 
hibits in response to any perturbing po tential can be computed 
(see, for example. iLeconte et al.ll201 ll). The Love number is a 
measure for the level of central condensation of an object; ^2 
decreases as the degree of central condensation increases. A ho¬ 
mogeneous incompressible ideal fluid body has ^2 = 3/2. For a 
star like CoRoT-29, the value of k 2 ca n be estimated from canon - 
ical stellar models (see, for example. IClaret & Gimene3ll995h . 
Considering the v sin / value and the mass and orbital distance of 
the companion, the expected value of ^2 for CoRoT-29 would 
produce a modest distortion, resulting in a J 2 that would be 3 
to 4 orders of magnitude smaller than our estimated value. In 
other words, the star is not rotating fast enough and the plan¬ 
etary companion is neither massive enough nor close enough to 
significantly distort the stellar shape. However, it is expected that 
the flattening may assume high values when the star is rotat- 
ing differentially, wi t h the angular velo city decreasing outwards 
(iJackson et al.ll20()^ IZahn et al.ll2010t) . In turn, in the assump¬ 
tion of shellular rotation, a stellar rotation regime where the an- 
gular velocity is c onstant on level surfaces but varying in depth, 
(IZahn et al.ll2010l) has shown that J 2 also increases compared 
to the case of uniform rotation. Assuming an internal differential 
rotation with a centre-to-surface ratio of 4 in angular velocity, the 
J 2 of a 7Mo rapidly rotating star like Archenar can exceed 0.01, 
which is about twice the value it assumes under the hypothesis 
of uniform rotation. It is not straightforward to infer how such 
a result could be extrapolated to a much less massive star and 
with slow surface rotation such as CoRoT-29. The star seems to 
be distorted and it is apparently a robust result, observed from 
space and independently confirmed from ground. The gravity 
darkening model provides the most accurate and reasonable ex¬ 
planation, although a dedicated modelling out of the scope of 
this paper would be required. 

5.2. Star-planet interaction 

In the previous section, we have mentioned that the stars seem 
to rotate faster than expected for their evolutionary stage. In the 
case of CoRoT-29, for example, if magnetic braking is acting in 
this KOV star, and it is expected to be, the rotation should be 
slower. So, either the system is very young, or we have to find 
a mechanism to speed up the stellar rotation One way to speed 
up rotation is tidal interaction between the star and the planet. 
But the planets CoRoT-28b and CoRoT-29b may not be very 


^ Note that ^2 is twice the apsidal motion constant, o ften indicated with 
the same symbol as in, e.g. lClaret & Gimened lll995h . 
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lag (days) 

Fig. 18. Autocorrelation function of the light curve of CoRoT-29. 



Fig. 19. Shape of the stellar surface of CoRoT-29 in the gravity darken¬ 
ing model (thick line) compared with a spherically symmetric star (thin 
line). 


efficient for that reason because of their orbital properties and, in 
the case of CoRoT-29b, for the high inclination of its orbit over 
the stellar equator ( for a review of star- planet interaction in the 
CoRoT context, see lPatzold et al.ll2012h . 

One hypothesis is that a second planet was present that ac¬ 
celerated the stellar rotation before falling into the star. Under 
this assumption, two planets formed simultaneously far away 
from the star (from the current position of the transiting planet). 
Planet-planet scattering could have moved then one planet in 
an inner orbit, very close to the star, and an outer planet in a 
larger, eccentric orbit. Tidal interactions would have made the 
inner planet fall into the star, accelerating its rotation rate and 
enriching its atmosphere with heavier elements such as lithium, 
while the outer planet would have subsequently circularized its 
orbit until the conditions that we see today. 

To assess the viability of this kind of a scenario, we investi¬ 
gated the behaviour of such a hypothetical planet in the neigh¬ 
bourhood of the star. Although it is impossible to find observa¬ 
tional evidence for this former planet, it is still useful to explore 
plausible mechanisms that can explain the current properties of 
CoRoT-29. 

To reduce the number of free parameters, we assumed that 
the planet was a twin of CoRoT-29b in a low-eccentricity orbit. 
Figure |20] shows the evolution of the orbit of such a planet and 
the rotation of the star under the joint action of the magnetic 
braking of the star and the tides on the star caused by the planet. 



TIME (Myr) 

Fig. 20. Evolution of the period of rotation of the star in days. The 
red curve is the orbital period of the hypothetical falling planet. Initial 
semi-major axis as indicated. Initial eccentricity: e=0.10. Initial rotation 
periods: 2, 5, and 12.8 days. 

Regardless of the initial rotational rate of the star, all simulations 
show that the stellar rotation quickly increases to 15 to 20 days. 
At the same time, the semi-major axis of the planet is decreasing. 
At some point the tidal torque of the planet on the star becomes 
important and this accelerates the stellar rotation rate. By the 
time the planet has collided with the star, the rotational period 
has been reduced to a few days. After the collision, the stellar 
rotation rate brakes to its present value. 

The models used in th e evolution simulation s are magnetic 
braking as determined by iBouvier et akl (1 19971) for stars with 
masses between 0.5 and 1.1 M© and a creep tide (iFerraz-Mellol 
12013l) due to a relaxation factor y = lOi '. This choice cor¬ 
responds to the mid of the tidal dissipation range det ermined 
from t he statistical st udy of stars hosting h ot Jupiters bv iHansenI 
(1201 Ol) (3-25 see lFerraz-Melloir2013l) . The great remaining 

difficulty is to discover how the existing planet CoRoT-29b could 
exist in same time as its falling sibling. The proximity of both 
orbits is prone to creating great instability; maybe the high incli¬ 
nation of CoRoT-29b is a consequence of this instability. 

6. Summary 

We have discovered and characterized two new planetary sys¬ 
tems: CoRoT-28b and 29b. 

CoRoT-28b belongs to the small population of hot-Jupiters 
orbiting evolved stars, which is interesting from the point of 
view of stellar and planetary evolution. It has a mass between 
that of Saturn an d Jupiter and does not seem to be inflated (see 
the discussion in lEnoch et alJl201^ . There are a couple of open 
questions about the characterization of CoRoT-28, in particular 
regarding its rotation. We have measured a higher than expected 
V sin i and its lithium content is higher than expected for the star’s 
evolutionary status. These open questions might reveal important 
information about the past of the system. For example, there is 
the question of whether there was a former planet that fell into 
the star and thus accelerated the stellar rotation and enriched 
its atmosphere with heavy elements. Unfortunately, we cannot 
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confirm this hypothesis because of the l ack of observational evi- 
dence. We used the method described in iParviainen et al.] (1201 3l) 
to search for planet occultations (secondary eclipses) in the light 
curves, but we did not hnd any signihcant signal. For CoRoT- 
28b we obtain a maximum star-planet flux ratio of 9.5% (99th 
percentile of the marginal posterior). This upper limit is not use¬ 
ful for constraining the planet’s albedo or brightness temperature 
in any useful way, since the highest physically feasible value for 
/ is ~ 2%. 

CoRoT-29b is a very interesting system. The strong asym¬ 
metry of the transit light curve is only comparable to transits of 
planets orbiting very young, fast rotating stars, but the star is 
neither young nor rotating fast. The extreme value of the de¬ 
formation, compared to stars of its mass range, is not under¬ 
stood and deserves more attention. We have not found at this 
point a satisfactory answer for the origin of the asymmetry in 
the transit shape, but we encourag e the study of thes e system 
with other instruments (CHEOPS, iBroeg et^ 1201 31 PLATO 
2.0. lRauer et al.llToi4 ground-based surveys) to conhrm or dis¬ 
prove CoRoT observations and solve the puzzle. 
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Fig. B.l. Flow chart for the gravity darkening modelling process of CoRoT-29b. 


Appendix A: Radial velocity data 

Appendix B: Flow chart for gravity darkening 
modelling 
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Table A.l. CoRoT-28 SOPHIE radial velocities, their errors, and bisec- Table A.4. CoRoT-29 HARPS radial velocities, their errors, and bisec¬ 
tor spans. tor spans. 


BJD 

RV 

[kms^'] 

Error 

[kms^'] 

Bis 

[kms^'] 

BJD 

RV 

[kms^'] 

Error 

[kms^'] 

Bis 

[kms^'] 

2455 801.34001 

-76.751 

0.023 

-0.055 

2456097.73381 

-64.131 

0.052 

-0.014 

2455 804.41446 

-76.849 

0.013 

-0.051 

2456098.70664 

-64.371 

0.046 

0.022 

2455 810.35554 

-76.782 

0.021 

0.012 

2456099.63447 

-64.199 

0.103 

0.030 

2455 811.32492 

-76.753 

0.013 

-0.019 

2456101.81698 

-64.303 

0.035 

-0.064 

2455 832.30855 

-76.723 

0.014 

-0.016 

2456102.80688 

-64.097 

0.056 

-0.177 

2456 063.56528 

-76.771 

0.029 

0.054 

2456115.66034 

-64.302 

0.077 

0.105 

2456 064.54259 

-76.820 

0.015 

-0.086 

2456116.71144 

-64.186 

0.026 

-0.076 

2456 071.54853 

-76.714 

0.017 

-0.049 

2456117.67788 

-64.083 

0.021 

-0.041 

2456 072.57780 

-76.718 

0.010 

-0.035 

2456119.61671 

-64.151 

0.051 

-0.091 

2456 100.49337 

-76.782 

0.012 

-0.128 

2456149.61905 

-64.255 

0.080 

-0.130 

2456103.48011 

-76.716 

0.010 

0.006 

2456150.63533 

-64.213 

0.037 

-0.095 

2456121.38285 

-76.810 

0.024 

-0.142 

2456151.63052 

-64.105 

0.045 

0.059 

2456 123.47526 

-76.671 

0.018 

-0.048 

2456154.57617 

-64.053 

0.038 

-0.082 

2456 124.46350 

-76.593 

0.019 

0.016 

2456159.55259 

-64.096 

0.056 

0.008 

2456 125.42640 

-76.671 

0.013 

-0.056 

2456161.53567 

-64.274 

0.037 

-0.087 

2456132.46267 

-76.844 

0.018 

0.010 

2456455.87546 

-64.189 

0.102 

-0.230 

2456133.40017 

-76.778 

0.011 

-0.000 

2456460.77835 

-64.325 

0.099 

-0.130 

2456134.43837 

-76.720 

0.012 

0.023 

2456508.68337 

-64.239 

0.041 

0.017 

2456 135.41925 

-76.709 

0.018 

-0.005 

2456511.54785 

-64.228 

0.031 

-0.130 

2456 149.44334 

-76.731 

0.017 

-0.037 

2456514.68026 

-64.315 

0.048 

-0.145 

2456 152.42582 

-76.797 

0.016 

-0.083 





2456153.41986 

-76.796 

0.010 

-0.089 





2456154.34156 

-76.750 

0.022 

-0.037 





2456 156.34663 

-76.732 

0.016 

-0.058 





2456 157.37367 

-76.765 

0.013 

-0.031 






Table A.2. CoRoT-28 HARPS radial velocities, their errors, and bisec¬ 
tor spans. 


BJD 

RV 

[kms-'] 

Error 

[kms-'] 

Bis 

[kms-‘] 

2456116.75525 

-76.7574 

0.0070 

-0.018 

2456118.75516 

-76.6534 

0.0076 

0.013 

2456152.54246 

-76.775 

0.013 

-0.034 

2456158.59125 

-76.7564 

0.0085 

-0.008 

2456 159.50759 

-76.7094 

0.0088 

-0.005 

2456161.49239 

-76.6614 

0.0083 

-0.057 

2456 514.54510 

-76.6742 

0.0085 

-0.028 

Table A.3. CoRoT-28 FIES radial velocities, their errors, and bisector 
spans. 

BJD 

RV 

[kms^'] 

Error 

[kms^'] 

Bis 

[kms^'] 

2456103.51866 

-76.800 

0.018 

0.004 

2456104.55555 

-76.862 

0.022 

0.001 

2456 105.56260 

-76.914 

0.017 

-0.033 

2456107.54147 

-76.844 

0.020 

-0.032 

2456117.47300 

-76.893 

0.011 

-0.008 

2456118.42105 

-76.810 

0.017 

-0.012 

2456119.50754 

-76.831 

0.015 

-0.009 

2456 120.55636 

-76.888 

0.015 

-0.034 

2456121.42907 

-76.897 

0.021 

-0.015 

2456 122.40610 

-76.876 

0.021 

-0.020 
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